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1. Introduction 

The study of an extended Higgs sector appears particularly adapted to the capabilities of a 
lepton collider of sufficient energy and luminosity. e + e~ and pL + pL~ collisions provide us with 
pair production of H + H and H A bosons. If photon collisions are available in a 77 collider 
setup, scans of the H° and A peaks can also be performed [p. In this paper we focus on the 
study of charged Higgs bosons at an e + e~ linear collider in supersymmetric scenarios. Since the 
heavy Higgs states generally decouple from the other Supersymmetric particles, most of, if not 
all, the results derived for the Minimal Supersymmetric extension of the Standard Model (MSSM) 
can be also applied to a non-Supersymmetric extended Higgs sectors, such as that of the two-Higgs 
doublet model (2HDM) or the Higgs triplet model [Q] . In the decoupling limit, with Ma > > Mz and 
the masses of the heavy Higgs particles much larger than those of the other states, the fundamental 
quantities to be determined in the charged Higgs study are the masses and the widths of these Higgs 
states and their decay branching fractions. The mass of one of the heavy Higgs bosons, generally 
taken to be the CP-odd A state, is a fundamental theory parameter. The A , H and ^ bosons 
are almost degenerate in mass in the decoupling limit and this needs to be verified experimentally. 
The heavy Higgs decay branching ratios are sensitive to another fundamental theory parameter, 
the ratio of the vacuum expectation values of the two Higgs doublets tan/3 = V2/V1. The total 
width of the heavy Higgs states is also important for determining tan/3. Finally, in SUSY models 
the contribution of sparticle loops may induce CP-violating asymmetries in the decay of charged 
Higgs bosons, which can be searched for at a linear collider. From an experimental point of view, 
the large jet multiplicity and event complexity of processes such as e + e~ — > H + H — > tbtb and 
W + h°W~h° make them excellent benchmarks for the detector granularity and its ability to perform 
accurate kinematic reconstruction and jet flavour tagging in an high particle density environment. 

2. Charged Higgs Production and Simulation Studies 

Charged Higgs production in e + e~ collisions proceeds mostly through s-channel pair produc- 
tion, e + e~ — > H + H . The production cross section at a given yfs energy depends on the charged 
Higgs mass and tan/3. Typical values range between ~20 fb for M H ± = 250 GeV at = 0.8 TeV 
and 0.5 fb for M H ± = 1140 GeV at -y/s = 3 TeV, due to the P-wave suppression near threshold (see 
Figure [I]). The single boson production e + e~ — > H~T + v r + c.c. []3| |J] and e + e~ — > H tb + c.c. 
[§] processes, which give access to charged Higgs production beyond the kinematic limit for pair 
production, have also been considered. Their cross sections scale as tan 2 /3. However, at ^(0.01 fb) 
they are too small to be exploitable in the interesting region M H ± > y/s/2. Production of charged 
Higgs bosons in t decays, t — > H + b has also been studied at a linear collider, using also the fact 
that the polarisation of the tau leptons from the subsequent H ± — > x ± v t decay is opposite to that of 
those originating from W ± — > t^v-i- for isolating the signal [^|. However, the Tevatron [^, |8]|, and 
soon the LHC, data are closing the region of interest for this process. 

In the mass range of interest for a linear collider, the main decay processes are H + — > tb, 
H + — > T + V T and H + — > W + h . The tbtb final state is dominant and it has been considered for 
most of the searches and mass measurement studies. Simulation studies of charged Higgs bosons 
have been performed at 0.5 TeV for M H ±= 140 and 180 GeV [g], at 0.8 TeV for M H ±= 200 and 
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Figure 1: Higgs pair production: cross section for e + e 
V5 = 0.5, 1.0ad3.0TeV. 



H + H at tan j3 = 20 as a function of M H ± for 



300 GeV [|10] |1_1|] and at 3.0 TeV for M H ±= 700 [jl2|, |13p, 900 qi4p and 1 140 GeV. The SM tbtb 
irreducible background is small, with a production cross section of 3.2 fb at 0.8 TeV and 1.2 fb 
at 3 TeV and a flat distribution in the tb invariant mass. The event reconstruction starts with the 
identification of two top quarks in events with large jet multiplicity and no significant missing 
energy. At energies below 1 TeV the boost of the top is small and the decay t — > W + b can be 
explicitly reconstructed. At multi-TeV energies and M H ± of 0.7 - 1.2 TeV, the top quark boost 
is such that it can be reconstructed as a single jet. Jet flavour tagging is essential since the four 
b jets are a distinctive signature of the WbbWbb final state. Once the four-parton final state is 
reconstructed, the t and b pair which minimises the mass difference between the two tb systems is 
selected. Given the large jet multiplicity, the four b tags and the need to reject leptonic W decays 
for an accurate t energy measurement, the reconstruction efficiency is of the order of just a few 
percent. 
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Figure 2: Mass reconstruction: tb invariant mass after kinematic fitting for signal events with M H ± 
200 GeV (left) and 300 GeV (right) at y/s = 0.8 TeV (from [M). 



The tb invariant mass resolution is 16 GeV for M H ± = 180 GeV at y/s = 0.5 TeV [g], 17 GeV 
for M H ± = 300 GeV at y/s = 0.S TeV [[Kj, 55 GeV for M H ± = 700 GeV Q, |13|] and 72 GeV for M H ± 
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= 1 140 GeV, both at y/s = 3.0 TeV. The mass resolution is due to the finite jet energy resolution, to 
missing energy from escaping neutrinos and to confusion in the association of individual particles to 
jets, due to the large multiplicities. This mass resolution can be improved by applying a kinematic 
fit to the reconstructed hadronic final states. Kinematic fits applied in simulation studies impose 
energy and transverse momentum p x =p y = 0,E±\p z \ = y/s, accounting for beamstrahlung photons 
radiated along the beam axis and an equal boson mass condition Mj\p, = Af/3;4- This improves 
the mass resolution by a factor of ~2 to 8.5 GeV for M H ± = 300 GeV at y/s = 0.8 TeV (see 
Figure @) [|T§], 25 GeV for M H ± = 700 GeV at y/s = 3.0 TeV |g ^ and 45 GeV for M H ± 
= 1 140 GeV. However, these resolutions are larger than the expected natural width at these 
masses. 



3. Discovery Reach and Present Predictions 

Since the SM background is small and flat at large values of the tb invariant mass, a significant 
signal of charged Higgs pair production can be established at an e + e~ linear collider, even with 
small cross sections and the sensitivity extends almost to the kinematical threshold for pair produc- 
tion by combining different decay final states, irrespective of the value of tan j8. At y/s = 3 TeV, the 
e + e~ — > H + H process is observable for charged Higgs boson masses up to 1.25 TeV (see Fig- 
ure ||) and exceeds the sensitivity of the process e + e~ — > H°A° which drops below 1 TeV at tan j8 ~ 
5 []12|]. This sensitivity can be compared to the expectations for the heavy Higgs mass values of 
new physics models within the constraints offered by previous searches, low energy data and, pos- 
sibly, astrophysical data, such as relic dark matter density. There have been extensive efforts to 
determine the regions of new physics model parameter space most likely to be realised in nature 
given these data. Most of these efforts concern some specific implementation of the supersymmet- 
ric extension of the Standard Model. The constraints most relevant here are those from the light 
Higgs boson mass limits, rare B decay rates and dark matter relic density. Given the large number 
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Figure 3: Likelihood marginalised to the - tan /3 plane with the contour for the anticipated sensitivity to 
H + H (dashed line) and H°A° (continuous line) pair production at a 3 TeV linear collider (plot modified 
from [p5[, sensitivity contours from [|T2|]) 
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of parameters of the MSSM, these studies are best performed in constrained versions, where the 
gaugino and sfermion masses are assumed to be unified at the GUT scale. This is generally known 
as the constrained MSSM (cMSSM). The physical masses are calculated using the renormalisation 
group equations (RGEs) from a set of five fundamental parameters: the universal scalar, mo, and 
gaugino, m-in, masses, tan/3, the universal trilinear coupling, Ao and the sign of the coupling /x. 
Different statistical techniques have been employed. A sampling of the full five-parameter space 



was performed in the study of ref. [ |15[ ] using a Markov chain Monte Carlo method [ ]16[ ] and taking 
into account limits from particle searches, £l x h 2 and low-energy constraints. Of particular interest 
here is the two-dimensional map of the Ma - tan/3 plane obtained, which is shown in Figure || 
where the sensitivity of a 3 TeV linear collider to H ± and H A pair production is superimposed. 
Variants of the cMSSM either add further constraints or relax them. In mSUGRA there are fixed 
relations between the bi-linear and trilinear couplings parameters and the scalar mass and the ad- 
ditional condition that the gravitino mass is equal to the scalar mass, m 3 / 2 = »iq, thus reducing the 



number of free parameters to just three [|7j|. In the non-universal Higgs mass model (NUHM) the 



supersymmetric contributions to the Higgs masses are allowed a different value [|18|]. In this model, 



predictions for the M^-tan /3 parameters have been obtained in the study of ref. [^9Q. A recent fre- 
quentist analysis, based on LEP-2 limits, electro-weak and flavour data, g — 2 and il^/j 2 results, 
considers all these models and determines likelihood functions for SUSY observables. The 95 % 
C.L. interval for the charged Higgs boson mass is found to be in the range 350 < M H ± < 750 GeV 
for the cMSSM and in the range 100 < M H ± < 600 GeV for the non-universal Higgs mass ex- 
tension [^]. A similar study in the mSUGRA model finds the interval 600 < M H ± < 1200 GeV, 
again at 95 % C.L. plj). These results show that present data place only a very loose constrain on 
the masses of the heavy Higgs particles and highlight the interest in pursuing detailed studies for 
charged Higgs bosons at the LHC and at future lepton colliders over a broad mass range up to, and 
possibly beyond, 1 TeV. 

3.1 Mass Determination 

The mass M H ± is determined by a fit to the two-parton invariant mass distribution after kine- 
matic fitting. The signal peak can be modelled by a Breit-Wigner function convoluted with a 
Gaussian resolution term. In general, the experimental resolution is larger than the H ± natural 



width. Alternatively, a template method, as that adopted in the study of ref. [12], can be used to 
build a % 2 for the mass fit. The SM background is essentially flat below this peak in the mass re- 
gion of interest. The residual H Q A° — > bbbb contribution gives a peaking background in the signal 
region, since the masses of the charged and heavy neutral Higgs bosons are expected to be almost 
degenerate. However, the e + e~ — > H + H cross section is larger than that for e + e~ — > H°A° and 
the t tagging reduces this background. The invariant mass distributions for signal and background 
at sfs of 0.8 TeV and 3.0 TeV are shown in Figure |] 

3.2 Decay Studies 



Modes other than the dominant tb decay have also been studied. In particular, the study of [ 12] 
has considered the leptonic decay H + — > T + V T in the mixed channel tbzv. Signal events can be 
selected using the standard reconstruction of the hadronic decay and the transverse mass of H + — > 
T + V T . The reconstruction efficiency, around 0.045, is comparable to that of the fully hadronic 
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Figure 4: Mass reconstruction: tb invariant mass for signal and background events reconstructed at -Js = 
0.8 TeV (left) (from [[]]]) and 3.0 TeV (right). 



Table 1: Summary of the results on the reconstruction efficiency and the determination of M H ± at various 
centre-of-mass energies. 
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145 


0.5 


0.5 


cs 


0.15 


0.006 






200 


0.8 


0.5 


tb 


0.02 


0.002 




300 


0.8 


0.5 


tb 


0.04 


0.004 
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702 


3.0 


3.0 


tb 


0.02 


0.007 
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12] 


880 


3.0 


3.0 


tb 


0.05 


0.008 


[ 


14] 


906 


3.0 


3.0 


tb 


0.07 


0.006 


Preliminary 


1136 


3.0 


3.0 


tb 


0.05 


0.007 


Preliminary 



mode. The determination of tan j3 will be essential to constrain phenomenology and to relate the 
extended Higgs sector to cosmology through dark matter. The precision of a linear collider is 
expected to provide essential inputs [^2]]. It is possible to constrain tan/3 by determining the H°, 
A , H + widths and decay branching ratios H°, A — > bb, t + t~ , H + — > tb, TV. In particular, the 
determination of the charged Higgs bosons decay yields to tb and cs hadronic and to tv t leptonic 
final states offers a good opportunity to determine tan /3 as long as its value is not too large [|3|, |2|]. 
At very large values of tan j8 complementary sensitivity is provided by the determination of the 
boson width as shown in Figure ||. 

In presence of CP violation, the scalar h° and H° bosons and the pseudo-scalar A are mixed, 
In this case the A is no longer a mass eigenstate and the H ± mass should be used instead of Ma for 
parametrising the MSSM [24]. CP violation in the new physics sector can manifest itself in decay 
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Determination of tan/S: Vi~=500 GeV, L=2000 fb 
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Figure 5: tanj3 determination with H ± : accuracy as a function of tanj3 forM H ± = 200 GeV (from [p3|]) 
rate asymmetries such as 

8 cp = T{H + ^ff')-nn-^ff) 
ff T (h+ -> //') + r(H- -> //') 

In particular, in MSSM CP asymmetries are expected to be mostly due to squark loops in H + — > tb 
decays p6| , p5| , ^]. The size of these asymmetries scale inversely with tan/3 from 0.20 to 0.02 
for 5 < tan j3 < 30. 8 CP can be determined using quark-anti quark tagging in hadronic top decays 
adopting a vertex charge algorithm for b and c jets and lepton charge in semileptonic t decays. 
The analysis is quite challenging due to the limited quark charge separation power and the small 
statistics. With typical cross sections of fb) and few ab _1 of integrated luminosity, a sensitivity 
to deviations of 8 CP from zero at the level of ~0.10 should be feasible. 

4. Conclusions 

A lepton collider of sufficiently high energy and luminosity can provide good accuracy in the 
determination of the mass (better than 1%), production cross section and decay branching fractions 
of charged Higgs bosons, through pair production virtually up to the kinematic limit. The accuracy 
afforded by an e + e~ linear collider will be essential for the interpretation of the nature of a non- 
minimal Higgs sector and understand its role in relation to relic dark matter. The study of the heavy 
Higgs sector is one of the important drivers towards high energy and high luminosity performance 
for a linear collider. At present the input from the LHC data is essential to define the machine and 
detector requirements. The ongoing activity in simulation studies for the ILC at 1 TeV and CLIC 
at 3 TeV will help clarifying requirements and physics potential based on realistic simulation and 
reconstruction including accelerator effects. 

Acknowledgements 

It is a pleasure to thank Tord Ekelof for the kind invitation to talk at this workshop and Arnaud 
Ferrari for useful discussion. 




7 



Charged Higgs Bosons at Linear Colliders 



Marco Battaglia 



References 

[1] J.I. Diana, LC-TH-2000-002. 

[2] K. Huitu, J. Laitininen, J. Maalampi and N. Romanenko, Nucl. Phys. B 598 (2001), 13. 
[3] S. Kanemura, S. Moretti and K. Odagiri, JHEP 0102 (2001) 011, [hep-ph/00 12030]. 
[4] S. Moretti, EPJdirect 4 (2002) 1, 

[5] B. Kniehl, M. Maniatis and M.W. Weber, arXiv: 1009.3929 [hep-ph]. 

[6] E. Boos, V. Bunichev, M. Carena, C.E.M. Wagner, ECONF C050318:0213 (2005) and 
arXiv:hep-ph/0507100. 

[7] T. Aaltonen et al [CDF Coll.], Phys. Rev. Lett. 103 (2009) 101803 [arXiv:0907.1269 [hep-ex]]. 
[8] V. Abazov et al [DO Coll.], Phys. Lett. B 682 (2009) 278 [arXiv:0908.181 1 [hep-ex]]. 
[9] A. Sopczak, CERN-PPE/93-197 
[10] M. Battaglia, A. Kiiskinen and P. Poyhonen, LC-PHSM-200 1-041. 

[11] R.D. Heuer, D. Miller, F. Richard, P. Zerwas (Editors), TESLA Technical Design Report, Part III 
(2001). 

[12] E. Coniavitis and A. Ferrari, Phys. Rev. D 75 (2007) 015004. 

[13] E. Coniavitis and A. Ferrari, PoS CHARGED2008 (2008) 022. 

[14] M. Battaglia, A. De Roeck, J. Ellis and D. Schulte (editors), CERN-2004-005. 

[15] B. C. Allanach and C. G. Lester, Phys. Rev. D 73 (2006) 015013 [arXiv:hep-ph/0507283]. 

[16] E.A. Baltz, P. Gondolo, JHEP 0410 (2004) 052 [arXiv:hep-ph/0407039]. 

[17] E. Cremmer, B. Julia, J. Scherk, P. van Nieuwenhuizen, S. Ferrara and L. Girardello, Phys. Lett. B 79 
(1978) 231. 

[18] H. Baer, A. Mustafayev, S. Profumo, A. Belyaev and X. Tata, JHEP 0507 (2005) 065 
[arXiv:hep-ph/0504001]. 

[19] J. R. Ellis, T. Hahn, S. Heinemeyer, K. A. Olive and G. Weiglein, JHEP 0701 (2007) 092 
[arXiv:0709.0098 [hep-ph]]. 

[20] O. Buchmuller et al. Eur. Phys. J. C 64 (2009) 391 [arXiv:0907.5568 [hep-ph]]. 

[21] O. Buchmuller et al. arXiv:101 1.61 18 [hep-ph]. 

[22] E. A. Baltz, M. Battaglia, M. E. Peskin and T. Wizansky, Phys. Rev. D 74, 103521 (2006) 
[arXiv:hep-ph/0602187]. 

[23] J. F. Gunion, T. Han, J. Jiang and A. Sopczak, Phys. Lett. B 565 (2003) 42 [arXiv:hep-ph/0212151]. 

[24] J. Ellis, J.S. Lee and A. Pilaftsis, Phys. Rev. D 76 (2007) [arXiv:0809.521 1 [hep-ph]]. 

[25] E. Christova, H. Eberl, W. Majerotto and S. Kraml, Nucl. Phys. B 639 (2002) 263 [Erratum-ibid. B 
647 (2002) 359] [arXiv:hep-ph/0205227]. 

[26] E. Christova, H. Eberl, W. Majerotto and S. Kraml, JHEP 0212 (2002) 021 [arXiv:hep-ph/021 1063]. 

[27] E. Christova, H. Eberl, E. Ginina and W. Majerotto, JHEP 0702 (2007) 075 [arXiv:hep-ph/0612088]. 



8 



